home *** CD-ROM | disk | FTP | other *** search
/ SGI Freeware 1999 August / SGI Freeware 1999 August.iso / dist / fw_plotutils.idb / usr / freeware / share / ode / henon.ode.z / henon.ode
Encoding:
Text File  |  1998-10-28  |  662 b   |  30 lines
  1. # You may run this example by doing:
  2. #
  3. #    ode < henon.ode | graph -T X -C
  4. # 
  5. # or alternatively, to get a real-time plot,
  6. #
  7. #    ode < henon.ode | graph -T X -C -x -0.6 0.6 -y -0.6 0.6
  8. #
  9.  
  10. # This example shows a simple center with complicated orbits:
  11. # the system investigated by Henon and Heiles.
  12. # The Poincare plot is a graph of p1(t) vs. q1(t).
  13.  
  14. # Initial condations:
  15. # stable:    p1 = 1/3, q1 = 1/4, p2 = 0.1293144, q2 = 1/5
  16. # unstable:  p1 = 0.1, q1 = 0.1, p2 = 0.467618,  q2 = 0.1
  17.  
  18. qone' = pone
  19. qtwo' = ptwo
  20. pone' = -qone - 2*qone*qtwo
  21. ptwo' = -qtwo - qone*qone + qtwo*qtwo
  22.  
  23. pone = 1/3
  24. qone = 1/4
  25. ptwo = 0.1293144
  26. qtwo = 1/5
  27.  
  28. print pone, qone
  29. step 0,250
  30.